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Abstract 

We reconsider the conflict between recent calculations of the semileptonic 
branching ratio of the B meson and the experimentally measured rate. Such 
calculations depend crucially on the application of "local duality" in nonlep- 
tonic decays, and we discuss the relation of this assumption to the weaker 
assumptions required to compute the semileptonic decay rate. We suggest 
that the discrepancy between theory and experiment might be due to the 
channel with two charm quarks in the final state, either because of a small 
value for m c or because of a failure of local duality. We examine the experi- 
mental consequences of such solutions for the charm multiplicity in B decays. 
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I. INTRODUCTION 



Because of the large energy which is released, the decay of a heavy quark is essentially 
a short distance process. This simple observation has led to much recent progress in the 
calculation of the inclusive decays of hadrons containing a heavy quark (l|-[7] • The method 
relies on the construction of a systematic expansion in the inverse of the energy release, 
given approximately by the heavy quark mass, and hence works most reliably in the bottom 
system. In fact, it is expected that certain features of inclusive bottom hadron decays may 
be reliably predicted with the accuracy of a few percent. 

Considerable attention has been paid to inclusive semileptonic and rare PJ^] B 

decays, both to total rates and to lepton and photon energy spectra. There is little con- 
troversy that these calculations rest on a firm theoretical foundation. However, it has been 
suggested to extend these methods to include nonleptonic decays as well [|,|| . This proposal 
has led to an intriguing conflict with experiment, as the predicted nonleptonic widths differ 
significantly from those which may be extracted from the measured semileptonic branching 
ratio of the B meson ||. In this calculation, the short-distance expansion has been carried 
out to third order in the inverse mass 1/mb, and a reasonable analysis leads the authors of 
Ref. || to the conclusion that it would be unnatural to find the source of the discrepancy 
in uncalculated terms of higher dimension or higher order in a s . 

It is the purpose of this article to reconsider this problem, in particular the assump- 
tions on which the computation is based. In Section II, we review the techniques used to 
treat inclusive decay rates, with an eye to emphasizing the differences between the theo- 
retical foundations underlying the calculations of semileptonic and nonleptonic decays. In 
Section III, we discuss the possible discrepancy between theory and experiment in the B 
semileptonic branching ratio. This might be resolved by an unusually small value for m c , or 
might involve the failure of the key assumption, "local duality" , underlying the calculation 
of the nonleptonic rate. In either case the enhancement of decays into final states with two 
charm quarks is a likely consequence. In Section IV we examine the implications of this for 
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the charm multiplicity in B decays, for which present data do not seem to support an en- 
hancement resulting from the b — ► ccs process. The unusual feature of the data on inclusive 
B decays is neither the semileptonic branching ratio alone, nor the charm multiplicity alone, 
but rather the combination of the two. Brief concluding remarks are given in Section V. 

II. THEORETICAL TECHNIQUES 

The weak decay of b quarks is mediated by operators of the form 

= J£J ei t, (2.1) 

where 

J^ = g l7 ^(l-7 5 )g 2 or - 7 >, (2.2) 

are fermion bilinears. The inclusive decay rate is given by a sum over all possible final states 
X with the correct quantum numbers, 

r~^(fi|C t |X)(X|C|5). (2.3) 

X 

In this article we adopt the notation that a generic B meson contains a b quark, rather than 
a b quark. The optical theorem may be used to rewrite Eq. ( [2.3|) as the imaginary part of a 
forward scattering amplitude, 

r~Im(£|T {£>*,£>} \B). (2.4) 

One then would like to use perturbative QCD to extract information about the time-ordered 
product appearing in Eq. ( |2.4j ). The extent to which this is possible is precisely the extent 
to which inclusive decay rates may be calculated reliably. 

In the case of semileptonic decays, one may follow a systematic procedure to justify the 
application of perturbative QCD JTJ. Up to negligible corrections of order diem and Gf, one 
may factorize the matrix element of the four-fermion operator, 
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{Xl{pz)v(p,)\J£J llx \B) = (X\ J£ \B)(£(p E )u(p p )\ J^O) , (2.5) 

and consider only the time-ordered product of the quark currents. One then finds an ex- 
pression in which the integral over the momenta of the leptons is explict, 

T~ J dydv qdq 2 L^v qj q 2 ,y)W^(v q,q 2 ) , (2.6) 

where L^ u is the lepton tensor and W^ v the hadron tensor. Here the momentum of the 
external b quark is written as p^ = m^iA The other independent kinematic variables are 
— Pe + Vp an d y = ZEi/rrib. It is convenient to scale all momenta by m&, so q = qjm^. 
The hadronic tensor is given by 

= Y.(B\j^\x)(x\r h \B) 

x 

= -21m (B\i J dxe*-T{ji?(x),Jjl(0)} \B) 

= -2ImT^. (2.7) 

One may perform the integrals in y, v-q and q 2 in Eq. (|2.6j ) to compute the total semileptonic 
decay rate, or leave some of them unintegrated to obtain various differential distributions. 

The doubly differential distribution dr / dy dq 2 is a useful case to consider. Here we must 
perform the integral over v ■ q, for y and q 2 fixed. The range of integration for v • q is given 
by (y + q 2 /y)/2 < v ■ q < (1 + q 2 — rh 2 )/2, where m q is the mass of the quark to which the 
b decays, and rh q = m q /mb. This integration is pictured in Fig. ||a, along with the analytic 
structure of T^ u in the v ■ q plane -Q ^ ne absence of a cut along the real axis in the 
region (1 + q 2 — rh 2 )/2 < v ■ q < ((2 + rh q ) 2 — q 2 — l)/2 is simple to understand in terms 
of the invariant mass pn of the intermediate hadronic state. Such a state may contain no b 
quarks, in which case it is subject to the restriction p 2 H = (rribV — q) 2 > m 2 (the left-hand 
cut), or it may contain bbq, in which case p 2 H = {m^v + q) 2 > (2mb + m q ) 2 (the right-hand 
cut). Except in the limit q 2 = g^ ax = (1 — m q ) 2 and m q — > 0, the two cuts do not pinch. 

*The discussion of the analytic structure of T^ u given in Ref. |fj] is erroneous. We thank B. Grin- 
stein and A.I. Vainshtein for discussions of this point. 
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In Fig. |T]a, we have already included only the imaginary part of by integrating over 
the top of the cut and then back underneath it. In general, along the physical cut will 
depend on v ■ q in a complicated nonperturbative way. We do not necessarily know how to 
compute in QCD in the physical region where there are threshold effects. However, we may 
use Cauchy's Theorem to deform the contour of integration until it lies away from the cut 
everywhere except at its endpoints, as illustrated in Fig. |l]b. Along the new contour, we are 
far from the physical region, and we may perform an operator product expansion for T^ v 
in perturbative QCD. Only far from any physical intermediate states is such a calculation 
necessarily valid. However, this is enough to allow us to compute reliably certain smooth 
integrals of T^ v by deforming the contour of integration into the unphysical region. That we 
can compute integrals of T^ u in perturbation theory in this way is the property of "global 
duality" . 

Unfortunately, the contour in Fig. |l] must still approach the physical cut near the end- 
points of the integration. This introduces an uncertainty into the calculation which cannot 
be avoided. Still, one has two arguments that this uncertainty is likely to be small. First, 
for large m&, the portion of the contour which is within Aqcd of the physical cut scales as 
Aqcd/^?) and thus makes a small contribution to the total integral. Second, if the energy 
release into the intermediate hadronic system is large compared to Aqcd, it is reasonable to 
expect that T^ u will be well approximated by perturbative QCD even in the physical region. 
This is because in this region the cut is dominated by multiparticle states, and hence the 
strength of the imaginary part of T^ u is a relatively smooth function of the energy. While 
new thresholds associated with the production of additional pions are found along the cut 
even in this region, their effect is small compared to the smooth background of states to 
which they are being added. 

This intuition, that for large enough energies one may perform the operator product 
expansion directly in the physical region, is "local duality" . While it is a reasonable property 
for QCD to have, it is obviously a stronger assumption than that of global duality. In 
particular, it cannot be justified by analytic continuation into the complex plane. Rather, it 
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rests on one's sense of how QCD ought to behave at high energies. It is clear, as well, that the 
energy at which local duality takes effect will depend on the operators which appear in the 
time-ordered product. Hence the fact that local duality appears to work at a given energy 
in one process, such as in electron-positron annihilation into hadrons, may be suggestive but 
does not prove that it should hold at the same energy in another process. 

To compute the inclusive semileptonic decay rate, then, one may use global duality except 
in a region along the contour of order AQCD/^b, where one must approach the physical cut. 
In this small region one must resort to local duality to justify the operator product expansion. 

Let us now turn to inclusive nonleptonic decays. Here there is no analogue of the factor- 
ization (|2.5|) which we had in the semileptonic case. Hence there is no "external" momentum 
q in which one may deform the contour away from the physical region, leaving one unable 
to use global duality in the transition to perturbative QCD. In this case, one is forced to 
invoke local duality from the outset if one is to argue that the the time-ordered product T^ v 
is computable. This clearly puts the calculation of inclusive nonleptonic B decays on a less 
secure theoretical foundation than that of inclusive semileptonic B decays. 

Nonetheless, we do not mean to assert that the assumption of local duality in nonleptonic 
decays is inherently unreasonable, merely that it is the least reliable aspect of the computa- 
tion. In fact, it is not entirely clear what it is reasonable to expect in this case. On the one 
hand, the energy released when a b quark decays is certainly large compared to Aqcd- On 
the other, the decay is initially into three strongly interacting particles (rather than into only 
one for semileptonic decays), and the energy per strongly interacting particle is not really 
so large. (Note that in the semileptonic case, the point at which the contour approaches the 
cut and local duality must be invoked is conveniently the point of maximum recoil of the 
final state quark, where local duality is expected to work best.) What we propose is that the 
comparison of the nonleptonic decay rate, as computed via the operator product expansion, 
with experiment be taken as a direct test of local duality in this process. As such, it is a 
probe of a property of QCD in an interesting kinematic region, and nonleptonic B decay 
well deserves the intense scrutiny which it has recently been accorded. 
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III. THE SEMILEPTONIC BRANCHING FRACTION OF B MESONS 



The experimental implications of inclusive nonleptonic decays of B mesons have recently 
been discussed in great detail by Bigi, Blok, Shifman and Vainshtein ||. Since the semilep- 
tonic branching ratio of the B is relatively well-measured, they use their calculation of the 
nonleptonic decay rate to predict this quantity. Their conclusion is that the semileptonic 
branching ratio which comes out of their computation is unacceptably high, corresponding 
to a nonleptonic width which is too low by at least 15-20%. In this section we will reconsider 
their analysis. 

The inclusive decay rate of the B meson may be divided into parts based on the flavor 
quantum numbers of the final state, 

r TOT = r(6 -> civ) + r(6 -> cud!) + r(& -> ccs') . (3.1) 

Here we neglect rare processes, such as those mediated by an underlying b —>■ u transition 
or penguin-induced decays. By d' and s' we mean the approximate flavor eigenstates (d! = 
d cos 6*i — s sin 9i, s' = d sin 9i + s cos 6\) which couple to u and c, respectively, and we ignore 
the effect of the strange quark mass. It is convenient to normalize the inclusive partial rates 
to the semielectronic rate, defining 

T(b -> cud') T(b -> ccs') 

Kud ~ 3 r( fo - ceu) ' cs ~ 3T(b ^ ceu) ' 1 j 

The full semileptonic width may be written in terms of the semielectronic width as 

r(6 -»• civ) = 3f(m T )T{b -»• ceu) , (3.3) 

where the factor 3/(m r ) accounts for the three flavors of lepton, with a phase space sup- 
pression which takes into account the r mass. Then, since the semileptonic branching ratio 
is given by Br{b — > civ) = T(b — > cIv)/Ttot, we may rewrite Eq. (|3.1|) in the form 

R ud + R cs = f{m T ) 1 ~ B ( [ ib ^ ■ (3-4) 

Br[b —> civ) 



The measured partial semileptonic branching fractions are [11,12 



Br(B -> XeP) = 10.7 ± 0.5% , 
Br(B -> X/iP) = 10.3 ± 0.5% , 

Br(.B -> Xri?) = 2.8 ± 0.6% , (3.5) 

leading to a total semileptonic branching fraction Br(b — > civ) of 23.8 ± 0.9%, with the 
experimental errors added in quadrature. Of the semileptonic rate, 11% comes from decays 
to r, corresponding to a phase space suppression factor f(rh T ) = 0.74, consistent with 
what one would expect in free quark decay |5|JI3|. If we substitute the measured branching 
fractions into the right-hand side of Eq. (|3.4|), we find 



R ud + R cs = 2.37 ±0.12. (3.6) 

We now compare this constraint with the theoretical calculations of R u d and R cs . 

The ratios R u d and R cs depend on the total rates T(b — > cez/), T(b — ► cud') and T(b — > 
ccs'). Each of these has a theoretical expansion in terms of a s (/i) and l/m&. Since corrections 
of order 1/m^ vanish and those of order l/mf are numerically expected to be at the few 
percent level |l|-4l|,||,|^Q] , we include here only the radiative corrections. Neglecting terms of 
order a 2 (/i), the expansions take the form 

r(6 - ceu) = r /(m c , 0) • jl - (vr 2 - ^ + 5 s . L (m c )) | , 

r(6 CMrf ') = r /(m c , o) • 3t7(//) |i - ^ - x + ^(^)) + j2 ^)} ' 

r(6 - ccs') = r o /(m c ,0) • 377(//)G(m c ) jl - (V - ^ + 5 cs (m c )) + J 2 (/i)| . (3.7) 

The prefactor To = G 2 F m\\ Vd,| 2 /1927r 3 will cancel in the ratios R u d and R cs , as will the 
charm quark phase space suppression J(m c ,0) |0|, to be discussed below. 

The radiative corrections have been computed analytically to order a s in the limit m c = 



0, and for semileptonic decays up to one numerical integration for general m c [[14]]. For 
semileptonic decays we absorb the correction due to ih c ^ into 5 s .i.(m c ) and present the 
numerical value of <5 s .i.(m c ) below. Finite charm mass effects for nonleptonic decays are 
absorbed into 8 u d{rh c ) and 5 cs {fh c ). Because these quantities have not been computed, we 
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present numerical results in the case of nonleptonic decays only for m c = 0. The expressions 
for F(b — > cud') and r(6 — > ccs') in Eq. ( |3.7| ) are somewhat more complicated than that for 
T(b — > ceP), due to renormalization group running between // = Mw and /i = m&. The 
leading logarithms are resummed into 77 = (2L^_ + L 2 _)/3, where |T5[ 



a. (Mi 



-i -6/23 



or, (Ml 



12/23 



(3.8) 



The subleading logarithms, which must be included if terms of order a a (fi) are also to be 
kept, are assembled into J2, 

J2 = ^rlT + 6l0g ^j2LlTL 2 : + 2 * 2LI + LI ' (3 ' 9) 



where p+ 



6473 



i, )(;(J(i — —0.51 and p_ = j|j|g = 1.47 arise from two-loop anomalous dimen- 
sions [16|. The factor 3 in Eq. (|3.7|) is for the sum over colors in the final state. Finally, 



there is an additional phase space suppression G{m c ) in T(b — > ccs') because of the masses 
of the two charm quarks. This factor is given by [13| 

I(rh r , rhr) 



G(m c ) 



l(m c ,0) 



(3.10) 



where 



l{x,0) = (1 -x 4 )(l -8x 2 + x 4 ) -24a; 4 log a: 



I(x, x) = Vl -4x 2 (l - 14x 2 - 2x 4 - 12x 6 ) + 24x 4 (l - x 4 ) log 



1 + Vl -4x r 



(3.11) 



1 - Vl - 4x 2 J 

In terms of the theoretical expressions (|3.7|) for the partial widths, the ratios take the form 



where 



Rud = P(ji) + $Pud(n, m c 



R cs = G(rh c 



P(n) + 5P cs (n,rh c 



(3.12) 



P(fi) = r](fi) 
5P ud (fi,m c ) = r)(n) 



1 H h J2W 
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2a g (/i) 
3tt 



4.1. (^ c ) - bud{rh c ) 



5P cs (fi,rh c ) = 77(/i) 2ns(/,) 



4.1. (m c ) - 5 cs (m c ) 



(3.13) 



3tt 

parametrize the radiative corrections. As emphasized in Ref. |1, if the theoretical expressions 



( |3.12j ) are inserted, then Eq. ( |3.6|) is not well satisfied. For example, if one simply takes the 
reasonable values /i = m b = 4.8 GeV, m c = 1.5 GeV, = 180 MeV and 5P u d = 5P CS = 0, 
then P(fi) = 1.27, G{fh c ) = 0.36 and the left-hand side of Eq. fl3.6| ) is only 1.73. We are 
thus tempted to push the uncertainties in the calculation as far as is reasonable, in order to 
see how much of the discrepancy can be resolved within the context of the operator product 
expansion. 

The largest uncertainty in the theoretical expression for R u d + R cs comes from the choice 
of the charm and bottom masses. Up to certain ambiguities which have recently been 



discussed [ 1T|| , within perturbation theory these masses should be taken to be the pole 
masses These masses have not been determined with much precision. However, 

within the heavy quark expansion, the difference between m c and m b is much more precisely 
known, in terms of the spin-averaged D meson and B meson masses: 

m b -m c = (M B ) ave . - (M D ) ave . = 3.34 GeV. (3.14) 

In what follows, we will hold mj — m c fixed, and consider variations of nib only. A reasonably 
conservative range for rrif, might be 4.4 GeV < m b < 5.0 GeV, which corresponds to 0.24 < 
m c < 0.33. In Fig. |2|, we plot G(rh c ) as a function of m bl using the constraint ( |3.14j ). In 
Fig. |3], we plot P{jj) for a variety of values of the QCD scale A^ |TT 



We start by considering R u d, for which the calculation is likely to be more reliable, since 
it is less sensitive to m c . There is uncertainty in the radiative correction P(n) from the 
choice of the renormalization scale fi. The usual choice fi = rnt is motivated by the fact 
that the total energy released in the decay is rrib- However, this energy has to be divided 
between three particles, so perhaps the appropriate scale is lower. For fi = 1.6 GeV ~ mb/3, 
a reasonable lower limit, and = 180 MeV, we find P{^) = 1.45, a modest enhancement 
over fi = 4.8 GeV. If Ag^ is taken as high as 220 MeV, we have P(/x) = 1.52, which 
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makes a small additional difference. The uncertainty in 5P u d is harder to estimate, since 
5 u d(rh c ) has not been calculated. However, one may extract 5 s ±(rh c ) by doing a numerical 



integration of the formulas in Ref. |TJ]]. For rh c = 0.30, we find 5 B ,i = —1.11, corresponding 
to (2a s (mb)/37r)5 s .i.(m c ) = —0.050. The magnitude of this correction grows approximately 
linearly with rh c , and for m c = 0.33, we have S s ± = —1.20. Hence the term is small and 
actually reduces Rud, although one might expect it to cancel in whole or in part against 
the term proportional to 8 U d(rh c ). What we can conclude at this point is that the error 
associated with ignoring the charm quark mass in the radiative corrections is likely to be no 
larger than ±0.05, and henceforth we will neglect this effect. 

The leading nonperturbative strong interaction corrections to R u d and R cs are char- 
acterized by the two dimensionless quantities K b = — {B(v)\ b v {iD) 2 b v \B{v)) /2m\ and 
Gb = {B{v)\b v g a G^ v a^ u b v \B{v)) /Arrvl. Because it breaks the heavy quark spin symme- 
try, the parameter Gb may be determined from the measured B* — B mass splitting, but 
the value of Kb is not known. Fortunately, Kb does not occur in the nonperturbative cor- 
rection to R u d- (Using the "smearing" technique of Ref. ||, this cancelation arises because 
F(b — > cud) and V(b — > ceu) have the same dependence on 777.5.) For R u d, then, we are more 
confident than for R cs that the nonperturbative QCD corrections are small. Note, however, 
that there is a contribution to the mass difference in Eq. (|3.14j) involving Kb and K c , which 
we have neglected. 

The above estimates lead us to the conclusion that with the effects we have included in 
the operator product expansion, it is difficult to avoid the upper bound R u d < 1.52. If this is 
true, then Eq. ( |3.6|) would imply R cs > 0.85. This can barely be achieved in the theoretical 
expressions we have given. If we vary 4.4 GeV < mj < 5.0 GeV, as suggested above, then 
0.27 < G(m c ) < 0.58. Estimating the radiative corrections as before, with A^L = 220 MeV, 
this suggests the upper limit R cs < 0.89, or R ud + R cs < 2.43. This is in agreement with 
experiment, but on the other hand, it requires us to push all the freedom in the calculation 
in the same direction, perhaps further than is reasonable. If one were to take the point of 
view that ji = 2.4 GeV ~ ra&/2 were the lowest reasonable value for fi, then one would have 
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the constraints R u d < 1.44, R cs < 0.83 and R u d + R cs < 2.27. If one were further to require 
rrib > 4.6 GeV, one would have R cs < 0.67 and R u d + R cs < 2.10. In this case, one might 
consider the discrepancy with experiment to be a more serious issue. 

Another possibility is that the relevant scale for the radiative corrections in the decay 
to two charm quarks is considerably lower than that for the final state with a single charm. 
Since the rest masses of the two charm quarks absorb approximately 60% of the energy 
available in the decay, the strongly interacting particles are not emitted with very large 
momenta. For example, the average energy of the strange quark in the decay b — > ccs, 
computed at tree level, is only about 1 GeV. With such a low energy the procedure of 
estimating the value of higher order QCD corrections by varying the subtraction point \i 
is of dubious value. In fact one might question whether any finite order of perturbation 
theory is adequate and whether threshold effects that cause a violation of local duality are 
important. 

It is evident from this discussion that nothing is particularly clear. Although the data on 
inclusive nonleptonic decays can almost be accounted for by squeezing the input parameters, 
one might feel a little nervous about the necessity of such a conspiracy. After all, as men- 

(5) 

tioned earlier the "reasonable" values fi = raj, = 4.8 GeV, m c = 1.5 GeV and A^ = 180 MeV 
lead to R u d = 1.27 and R cs = 0.46, far short of the mark. An enhancement of approximately 
40% in the nonleptonic rate is called for. If one were to require this effect to be found 
entirely in R cs , it would amount to more than a factor of two. While we are less inclined 
than the authors of Ref. to insist that something is amiss, it is nonetheless intriguing to 
consider the possibility that the data indicate an enhancement of the nonleptonic rate over 
and above what we have included in the operator product expansion. Where might such an 
enhancement come from? 

The simplest explanation would be that due to a failure of local duality, the inclusive 
nonleptonic decay rate is simply not calculable to better than 40% or so. This is certainly 
a discouraging explanation, in that if it were true then there would be very little one could 
say in detail about why local duality, and hence the calculation, had failed. One was simply 
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unlucky. On the other hand, this explanation may well be correct. While we expect local 
duality to hold in the asymptotic limit of infinte b quark mass, we have little to guide us 
in estimating how heavy the b quark actually needs to be in practical terms. In particular, 
it is not relevant to consider, at low orders in QCD perturbation theory, the size of a 
few subleading terms which appear in the operator product expansion itself. The matrix 
elements which appear in this expansion are sensitive to details of the B meson bound state, 
but they are explicitly not sensitive to resonance effects in the final hadronic state. 

If local duality fails, it could well fail differently in the T(b — > cud!) and Y(b — > ccs') 
channels. In fact, we would expect it to fail worse in the channel with two charm quarks, 
since we expect the final states to be characterized by lower particle multiplicity and be 
closer to the resonance-dominated regime. Local duality, by contrast, is applicable only 
in the regime where the effect of individual resonance thresholds is small compared to the 
almost smooth "continuum" of multiparticle states. On the other hand, the phase space 
suppression from the two final state charm quarks means that unless m c is unusually small, 
only thirty percent or so of the inclusive nonleptonic rate comes from the T(b — > ccs') 
channel. Hence, to account for an enhancement of the full nonleptonic rate by forty percent 
purely from b — > ccs' would require a dramatic failure of local duality in this channel. 

IV. EXPERIMENTAL CONSEQUENCES OF AN ENHANCEMENT OF R cs 

Either through a failure of local duality, or from an unusually small value for m c , or 
because of a combination of these effects, the value of R cs is likely to be near unity in order 
to account for the measured B semileptonic branching ratio. This corresponds to about 
one-third of B decays arising from the b — > ccs' process. One consequence of this is a large 
number of charmed quarks per B decay, 



We remind the reader that we have adopted the notation that a generic B meson contains 
a b quark, rather than a b quark. Using Br(B — > X c lv) = 23.8% and f(m T ) = 0.74 in 



n c — 1 + R, 



■cs 




(4.1) 
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Eq. fly) yields 

n c = 1.00 + 0.32 R cs , (4.2) 

which for the values of R cs necessary to explain the semileptonic branching ratio would 
indicate n c ~ 1.3. 

There are contributions to the experimental value of n c from charmed mesons, charmed 
baryons, and cc resonances. The number of charged and neutral D mesons per decay, 



summed over B and B, has been measured to be |19 



n D ± = 0.246 ±0.031 ±0.025, 
n D0 - f , = 0.567 ±0.040 ±0.023. (4.3) 

The branching ratio to Df mesons has not yet been determined, because no absolute D s 



(4.4) 



branching ratio has been measured. However, it is known that p0 | 

3 7% 

n n ± = (0.1224 ±0.0051 ±0.0089) — — — 

Ua [Br(D s — > (f)n) 

and the branching ratio for D s — > <f)ir is expected to be about 3.7%. 

We must include in n c twice the inclusive branching ratio to all cc resonances which 
are below DD threshold. The measured inclusive branching ratio to ip is (1-H ± 0.08)%, 
including feed-down from ip' and \ c decays [jnj. It is also known that Br(B — > ip'X) = 
(0.32 ± 0.05)%, Br{B -> X ciX) = (0.66 ± 0.20)% and Br{B -> r) c X) < 1%. Hence we 
expect that the inclusive B branching ratio to charmonium states below DD threshold is 
about 2%. 



The inclusive B decay rate to baryons is about 6% [0. While it is commonly believed 
that these baryons arise predominantly from the b cud' process, giving A C X final states, 
we argue elsewhere |21| that a large fraction of B decays to baryons actually arise from the 
b — ► ccs 1 process, which gives final states with both a charm baryon and an anticharm baryon, 
such as S C A C X. Evidence for this interpretation comes from the experimental distribution of 
A c momenta, which shows that the A c 's produced in B or B decay are recoiling against a state 

14 



with a mass greater than or equal to the mass of the H c [p^-p4|. This novel interpretation 



of B decays to baryons can be consistent with the measured A£ ± correlations if Br(S c — ► 
AX)/Br(A c -> AX) is large 0. 

Even if D decay to baryons predominantly gives final states with both a charm and an 
anticharm baryon, the data summarized above do not provide supporting evidence for a 
value of n c around 1.3. Given the uncertainties, however, such a large value for the number 
of charmed hadrons per B decay is perhaps not excluded. From our perspective the curious 
feature of the data on inclusive B decay is not the measured semileptonic branching ratio 
alone, but rather the combination of it with the data on charm multiplicity in these decays. 

In this paper we have neglected B decays that do not arise from an underlying b — > c 
transition. Other possible processes include the b — » u transition and contributions from 
penguin-type diagrams. While it is very unlikely that such sources contribute significantly 
to the nonleptonic decay rate, this assumption can be tested experimentally, if enough 
branching ratios can be measured. The fraction of B decays arising from the b — > c transition 
is given by the sum of the B branching ratio to charmonium states below DD threshold, the 
branching ratio to states containing at least one charmed baryon, and the branching ratios 
to the ground state charmed mesons, BriB — > D°X), Br(B — ■> D + X) and Br(B — > DfX). 
Note that the inclusive charm yields reported in Eqs. ( j4.3|) and ( j4.4|) are actually sums of 
branching ratios (for example, neglecting CP violation, n D ± = BriB — > D + X) + BriB — ► 
D~X)). However, it should be possible with enough data to extract the individual branching 
ratios themselves. 

For example, one could count the number of DD (or D D) events per BB event at the 
T(4£), 

n(DD) = (1 - b) Br(B -»■ DX) Br(B -> DX) 

+- [Br{B -> DX) Br(B -> DX) + Br(B -> DX) Br(B -> DX) 

— b 

w (1 - b)Br(B -> M) Br(B -> DX) + - Br(B -> DX) Br(5 -> DX) , (4.5) 

and 
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n{DD+) = {l-b) \Br(B -> DX) Br(B -> D+X) + £r(£ -> D+X) Dr(DT -> DX) 
+^ [Sr(S -> -DX) £r(D -> D+X) + Dr(D -> DX) Dr(D -> D+X) 
« (1 - 6) Dr(D -»• DX) Br(B -»• D+X) . (4.6) 



where 6 ~ 0.076 is the B — B mixing parameter pq| . It is defined as the fraction of B meson 
events which are BB or B B, and is measured directly from lepton-lepton sign correlations, 



N B b + AW + N e 



(4.7) 



N BB + n bb + N l+l - + N i+i+ + N t -t- 

Combining Eqs. ( |4.3| ) and ( [4.5[ ), we may extract BriB — > DX) and Br(B DX) sepa- 
rately, if we neglect CP violation and impose the constraint BriB — > DX) = Br(B — ► DX). 
Analogously, we may extract Br(B — > D+X) and BriB — > D~X). Another method for 
determining individual branching ratios would involve tagging the flavor of the D which 
produced the charmed hadron by measuring the charge of a hard primary lepton from the 
other D in the event. 

Invoking a large violation of local duality has some implications for the pattern of B 
meson decays which may be different from what would be expected if local duality held 
and an unusually small value of m c were used to explain the measured B semileptonic 
branching ratio. For example, a violation of local duality in the b — ► ccs channel could lead 
to quite different lifetimes for the B, B s and A&, differences which are small in the operator 
product expansion because they arise only from higher dimension operators. However, since 
the effective Hamiltonian for this process has isospin zero, the equality of the B° and B~ 
lifetimes would not be disturbed. Similarly, violations of local duality in b — > cud' could lead 
to unequal B° and B~ lifetimes. B decay event shapes can also provide a test of the free 
quark decay picture for the b — ► cud' decay channel |26| . 



V. CONCLUDING REMARKS 



We have examined whether the measured B meson semileptonic branching ratio can be 
explained within the conventional application of the operator product expansion, in which 
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operators of low dimension are kept and perturbative corrections are included to a few orders 
in a s . We have found that this scenario would require an unusually small value for m c . If 
instead the explanation lies outside the conventional application of the operator product 
expansion, then a failure of local duality in the b — > ccs' channel is the likely explanation 
for the discrepancy with experiment. In either case, we expect the number of charmed 
hadrons per B decay to be approximately 1.3. Unfortunately, the present data on charm 
multiplicities do not support such a large value of n c . From our perspective, the unusual 
feature of inclusive B decay is not the semileptonic branching ratio alone, nor the charm 
multiplicity alone, but rather the combination of the two. Together, they would seem to 
suggest a significant violation of local duality in the b — > cud' nonleptonic decay process. 
From a theoretical point of view, however, such a resolution would be somewhat unsettling, 
as it would indicate a breakdown in the computation of the nonleptonic decay rate in the 
region where it is expected to be the most reliable; we understand why such a conclusion 
was resisted by the authors of Ref. ||. Still, it remains an open possibility, indicating 
perhaps that the invocation of local duality in quark decay requires a considerably larger 
energy release than has been naively hoped or expected. Given the apparent difficulties in 
performing a reliable computation of the nonleptonic decay rate, then, the CKM matrix 
element V c b should be extracted from the B semileptonic decay width rather than from the 
B lifetime. The uncertainties in such an extraction arise primarily from the choice of mj 
and subtraction point and are discussed in detail in Refs. p7p8| . 
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FIGURES 

FIG. 1. Contours in the complex v ■ q plane, for fixed q 2 and y. The gap between the cuts 
extends for (1 + q 2 — rh 2 )/2 < v • q < ((2 + rh q ) 2 — q 2 — l)/2. The endpoints of the contour integral 
are at v ■ q = (y + q 2 /y)/2 ± ie. 

FIG. 2. The phase space suppression factor G(m c ), as an implicit function of m& with 
m b -m c = 3.34 GeV held fixed. 

FIG. 3. The radiative correction P{n). The upper curve corresponds to = 220 MeV, the 
middle curve to A^i = 180 MeV, and the lower curve to = 140 MeV. We take m b = 4.8 GeV. 
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Figure 1a 




Figure 2 
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